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ABSTRACT 

The X-ray transient XTE J1701-462 was the first source observed to evolve through all known subclasses of 
low-magnetic-field neutron star low-mass X-ray binaries (NS-LMXBs), as a result of large changes in its mass 
accretion rate. To investigate to what extent similar evolution is seen in other NS-LMXBs we have performed 
a detailed study of the color-color and hardness-intensity diagrams (CDs and HIDs) of Cyg X-2, CirX-1, 
and GX 13+1—three luminous X-ray binaries, containing weakly magnetized neutron stars, known to exhibit 
strong secular changes in their CD/HID tracks. Using the full set of Rossi X-ray Timing Explorer Proportional 
Counter Array data collected for the sources over the 16 year duration of the mission, we show that Cyg X-2 
and Cir X-l display CD/HID evolution with close similarities to XTE J1701-462. Although GX 13+1 shows 
behavior that is in some ways unique, it also exhibits similarities to XTE J1701-462, and we conclude that its 
overall CD/HID properties strongly indicate that it should be classified as a Z source, rather than as an atoll 
source. We conjecture that the secular evolution of Cyg X-2, Cir X-l, and GX 13+1—illustrated by sequences 
of CD/HID tracks we construct—arises from changes in the mass accretion rate. Our results strengthen previous 
suggestions that within single sources Cyg-like Z source behavior takes place at higher luminosities and mass 
accretion rates than Sco-like Z behavior, and lend support to the notion that the mass accretion rate is the 
primary physical parameter distinguishing the various NS-LMXB subclasses. 

Key words: accretion, accretion disks - stars: neutron - X-rays: binaries - X-rays: individual (CirX-1, 

Cyg X-2, GX 13+1, XTE J1701-462) 


1. INTRODUCTION 

In neutron star low-mass X-ray binaries (NS-LMXBs), a 
neutron star accretes matter from a low-mass companion star. 
The vast majority of NS-LMXBs contains weakly magnetized 
neutron stars, and these systems are traditionally divided into 
two main subclasses based on their correlated spectral and 
timing properties: the Z sources and atoll sources (Hasinger 
|& van der Klis[[1989| ). The names derive from the shapes of 
the tracks the sources trace out in the X-ray color-color and 
hardness-intensity diagrams (CDs and HIDs). The Z sources 
are the more luminous, with near-Eddington X-ray luminosi- 
tities, whereas the atoll sources are thought to in general have 
Lx < 0.5L Edd . Further divisions of these two main NS-LMXB 
types have also been made; in particular, the Z sources have 
been divided into the Cyg-like and Sco-like subtypes ( jKuulk-| 
|ers et al.|1997| and references therein). 

In addition to differences in X-ray luminosity, rapid- 
variability characteristics, and spectral properties (see, e.g., 
|van der Klis|20061 and references therein), the NS-LMXB sub- 
classes also differ with respect to behavior in the radio band 
(Migliari & Fender 2006|) a nd the rat es an d properties of type I 
X-ray bursts exhibited ( |Galloway et al.|2008J ). Understanding 
what physical factors underlie tins" variety in characteristics 
has been a long-standing problem in the study of X-ray bina¬ 
ries. It has long been clear that the mass accretion rate must 
play an important role (e.g., it seems evident that Z sources 
must in general accrete at a higher rate than atoll sources, 
given the former’s significantly higher luminosity), but ex¬ 
actly what that role is, and where other physical parameters 
enter the picture, has been debated. 

Another long-standing question concerns the physical na¬ 


ture of the motion of Z sources in the CD/HID. The Z sources 
trace out characteristic tracks in the CD and HID (consist¬ 
ing of the so-called horizontal, normal, and flaring branches) 
on typical timescales of hours to a few days. Initially, it 
was assumed that motion along the Z track was driven by 
changes in the mass accretion rate, with M increasing mono- 
tonically from the tip of the horizontal branch to the tip of 
the flaring branch (e.g., Hasinger et al.||l990] ), but later ob¬ 
servations cast doub ts on this, and other scenarios have been 
proposed (see, e.g., |van der Klis||2001[ |Homan et al.||2002, 
2007} |Church et al.| 2012| ). In addition to motion along the 


tracks these sources show another type of motion, usually on 
longer timescales (days to weeks), where the tracks can shift, 
primarily in the HID. This has been referred to as secular mo¬ 
tion/shifts/changes/variations. In most cases these shifts are 
small and do not lead to significant changes in the shapes of 
the tracks. However, a few sources have been known to ex¬ 
hibit much stronger secular changes where the tracks shift and 
change shape radically in both the CD and HID. 

A breakthrough in our understanding of the Z/atoll phe¬ 
nomenology and the nature of secular evolution came with the 
transient NS-LMXB XTE J1701-462, which over the course 
of its 2006-2007 outburst evolved through all subclasses of 
low-magnetic-field NS-LMXBs—from a Cyg-like Z source at 
the highest luminosities to a Sco-like one, followed by a phase 
in the atoll source soft state (during which type I X-ray bursts 
were seen), and ending with a transition to the atoll source 


hard state before returning to quiescence (|Lin et al. 2009b; 
Homan et al.|2010| ). |Lin et al.| ( |2009b| ) performed a detailed 


spectral analysis of the entire outburst, using data from the 
Rossi X-ray Timing Explorer ( RXTE ), and their results indi- 
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cate that changes in the mass accretion rate_were responsi¬ 
ble for the evolution of the source. Homan et al.| ( |2010| )— 


Table 1 

Total Exposure Times for Analyzed Sources 


behavior of the source implies that differences in the mass 
accretion rate can alone explain the existence of the various 
NS-LMXB subclasses, and that it is not necessary to invoke 

Source 

Exp. Time (Ms) 

XTE J1701-462 

2.71 

differences in other parameters, such as the magnetic field of 

Cyg X-2 

2.28 

the neutron star (Hasinger & van der Klis 1989; Psaltis et al. 

Cir X-l 

2.57 

1995) or the inclination angle of the system (Kuulkers & van 

GX 13+1 

0.58 

der Kiisjtqq^). 


(if any) evolution similar to that of XTE J1701-462 is ob¬ 
served in other sources. Comparisons of XTE J1701-462’s 
evolution at the low-luminosity end with atol l sources have 
been made i n |Lin et al.] ( |2009b| ), |H10[ and |Munoz-Darias| 


et al. ( [20141 ), showing that there the source behaved simi¬ 
lar to otherlow-luminosity (atoll) NS-LMXBs; in this paper 
we focus mainly on the high-luminosity (Z source) portion of 
XTE J1701-462’s evolution. To this end we performed a de¬ 
tailed analysis of the CDs and HIDs of three luminous neutron 
star X-ray binaries—Cyg X-2, Cir X-l, and GX 13+1—using 
all available RXTE Proportional Counter Array (PCA) data. 
Similar to XTE J1701-462, these sources are known to have 


shown strong secular evolution (e.g., Oosterbroek et al. 1995 

Kuulkers et al. 1996; Wijnands et al. 1997 

Shirey et al. 1999a 

Schnerr et al. 2003). Since we are mainly i 

nterested in sources 

include in our study 

that cross subclass boundaries, we do not 


persistent Z sources that show only mild secular shifts and 
no significant changes in the shapes of their tracks. Although 
small subsets of the RXTE data sets for Cyg X-2, Cir X-l, 
and GX 13+1 have been studied in several papers—usually 
with an emphasis on spectroscopy and/or timing analysis— 
we are not aware of any papers presenting a comprehensive 
study of the secular evolution of any of these sources using 
a large am ount of RXTE data, with the exception of |Schn-| 
| err et al.| (|2003|), who studied a large set of observations of 
GX 13+1 made in 1998 (see discussion in Section [T4|). How¬ 


ever, since then the amount of GX 13+1 data in the RXTE 
archive has grown considerably. We also note that Shirey 
and collaborators studied CDs and HIDs of Cir X-l using a 
numbe r of observations made in 1996 a nd 1997 in a series of 
pap ers ([Sh irey et a l.|1996|fl998[|[1999a[ see discussion in Sec¬ 
tion |3.3| )7T>ut these observations represent only a small frac¬ 
tion of the currently available RXTE data for this source. In 
addition to our analysis of Cyg X-2, Cir X-l, and GX 13+1 
we also perform some reanalysis of XTE J1701-462 in order 
to ensure complete consistency with the analysis of the other 
three sources and to facilitate comparisons between them and 
XTE J1701-462. We note that although Cir X-l may be a 
rare low-magnetic-field neutro n star high-mass X-ray binary 
(as further discussed in Section [33) we will for simplicity in 
general refer to all four sources studied in this paper as NS- 
LMXBs. 

The structure of the paper is as follows. In Section [2] we 
describe the general data analysis steps common to allfour 
sources, but leave the descriptio n of analysis specific to each 
source to Section [3] In Section |3.1| we briefly describe our 
analysis and present our results for XTE J1701-462. This 
analysis largely follows that previously performed by H10 
(and includes creating a sequence of CD/HID tracks that 
shows the secular evolution of the source); we therefore 
mostly refer to that paper for details. In this section we also 
tie this source’s behavior to the various source states/branches 
generally seen in the CDs and HIDs of NS-LMXBs. In Sec¬ 


tions T2f[3.4| we describe our analysis and present our results 
for Cyg X-2, Cir X-l, and GX 13+1, respectively. For each 
source we construct a sequence of CD/HID tracks (analogous 
to the one for XTE J1701-462) that illustrates its secular evo¬ 
lution; we also give a brief background on each of these three 
sources. In Section [4] we discuss our results, and in Section [5] 
we give a summary of our results and conclusions. 

2. DATA ANALYSIS 
2.1. Data Extraction 

We used the PCA data from all pointed RXTE observa¬ 
tions of XTE J1701-462, Cyg X-2, Cir X-l, and GX 13+1 
obtained during the lifetime of the mission. In Table [T] we 
list the total amount of useful data (i.e., data remaining af¬ 
ter the removal of X-ray bursts and bad data of any kind) for 
each of these sources. In the case of Cir X-l, however, the 
quoted exposure tim e includes data affected by local absorp¬ 
tion; see Section |33| Because of the large number of observa¬ 
tions (over 2300 individual ObsIDs), most of the data analysis 
steps described below were automated. 

The data were analyzed using HE AS OFT, versions 6.9- 
6.12, as well as locally developed software. We used the 
Standard-2 data, extracting counts in 16 s time bins, com¬ 
bined from all three xenon layers, from all active Propor¬ 
tional Counter Units (PCUs) at any given time. The data were 
filtered following the recommendations of the RXTE Guest 
Observer Facility (GOF)jjthis included the removal of data 
around PCU voltage breakdown events (but only from the rel¬ 
evant PCU in each case). We also corrected for dead time us¬ 
ing the standard procedure recommended by the RXTE GOF 
and subtracted background using the faint or bright back¬ 
ground model as appropriate in each observation, based on 
the average count rate during the observation after exclusion 
of any type I X-ray bursts. 

2.2. Burst and Particle Flare Removal 

We removed all data obtained during type I X-ray bursts; to 
identify bursts in an aut omated fashion we used a method sim¬ 
ilar to the one used by [Remillard et al.[ ( [2006| ). We note that 
although some bursts observed from Cyg X-2 and Cir X-l do 
not show cooling along their tails—as is usually observed in 
type I X-ray bursts—the origin of these faint burs ts is nev¬ 
ertheless very likely thermonuclear ( [Linares et al.|20lT] ) and 
they were removed. 

PCUs 0 and 1 lost their propane layers in 2000 and 2006, 
respectively, diminishing their ability to reject events due to 
particle background. Based on a comparison with data from 
PCU 2—which did not suffer from such particle flares and 
was nearly always on—these events were identified and sub¬ 
sequently removed from the data. We note that since PCU 2 

4 See http://heasarc.gsfc.nasa.gov/docs/xte/xhp_proc_ 
analysis.html 
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data exist for all times during which PCU 0/1 data were re¬ 
moved, we did not exclude data from any given point in time 
entirely, but simply reduced the number of PCUs contributing 
to a particular time bin in these cases. 


2.3. Response Correction 

From the Standard-2 data we extracted count rates in sev¬ 
eral different energy bands and used those to calculate colors 
for CDs and HIDs. The response of the RXTE PCUs evolved 
over the duration of the mission due to several factors, and 
these changes in response must be corrected for before data 
across the lifetime of the mission can be combined into a sin¬ 
gle CD or HID for a given source. For this correction we 
used observations of the Crab Nebula, which was observed 
on average ~3 times per month throughout the mission. For 
each Crab observation we extracted the average count rate in 
each of the energy bands of interest for each active PCU, us¬ 
ing the Standard-2 data. We then fitted the mission-long light 
curve (where each data point represented the average count 
rate during a single observation) for each PCU in each en¬ 
ergy band with a piecewise-linear function. We (somewhat 
arbitrarily) chose MJD 50800 (1997 December 18) as our ref¬ 
erence epoch and normalized all data points to that date, first 
for each PCU individually, and then in each case the count 
rates from PCUs 0, 1,3, and 4 to that of PCU 2. 

Recently it has become apparent that X-ray emission 
from the Crab Nebula shows significant intrinsic variability. 
|Wilson-Hodge et al.| (|201 1J) show that from 2001 to 2010 
the RXTE PCA count rate from the Crab (after correcting for 
changes in response) varied rather irregularly by several per¬ 
cent. The (response-corrected) PCU 2 variability is ^5% in 
the 2-15 keV band and ~8% in the 15-50 keV band. This 
can compromise our correction for the variation in the PCA 
respons e. Comparing our mission- long Crab light curves to 
those of |Wilson-Hodge et al.| ( |201 1| we can estimate the mag¬ 
nitude of the effect on our derived rates and colors. The rel¬ 
ative amplitude of the Crab count rate fluctuations gradually 
increases with energy, but otherwise the variability behaves 
in more or less the same way in the different energy bands. 
Due to this energy dependence of the variability it should have 
some effect on the colors. However, given the relatively small 
difference in the strength of the variability between adjacent 
energy bands, we expect that this variability will largely can¬ 
cel out in the colors. We estimate that shifts in our colors due 
to the Crab variability are at most ^l%-2% and in most cases 
significantly less; we therefore expect the influence of this on 
our results to be negligible (cf. dis cussi on about uncertainty 
due to counting statistics in Section [2~4] ). The intensity is af¬ 
fected more strongly; we expect that shifts there can possibly 
be as high as ^6% in the worst case. However, they are gen¬ 
erally much smaller—probably less than ~4% in almost all 
cases and typically in the 0%-3% range—and we expect the 
effect of this on our results to in general be negligible as well. 


2.4. Construction of CDs and HIDs 

For the creation of CDs (i.e., hard color versus soft color) 
and HIDs (hard color versus intensity) we used color defini¬ 
tions similar to those used in lHlQI we defined our soft color 
as the net counts in the 4.0-7.3 keV band divided by those 
in the 2.4-4.0 keV band, and our hard color as net counts in 
the 9.8-18.2 keV band divided by those in the 7.3-9.8 keV 
band. The intensity we used for the HID was the net count 
rate per PCU in the 2-60 keV band. Before creating our CDs 


1.4 


1.2 
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Figure 1 . CD and HID representing the entire RXTE PCA data set of 
XTE J1701-462. 

and HIDs we combined the counts from all active PCUs for 
each 16 s time bin (after performing the corrections and fil¬ 
tering described above) and then rebinned the data in a given 
observation in order to maintain a more uniform size in the 
error bars across different values in count rate. We set a min¬ 
imum of 16,000 counts in the 2-60 keV band (after back¬ 
ground subtraction) for each rebinned data point. The rebin¬ 
ning was done in an adaptive/dynamic fashion—i.e., we in 
general did not use a single value for the time binning factor 
over an entire observation, but instead allowed the factor to 
vary (while imposing our counts limit) to adapt to a possibly 
varying count rate. In a few cases we applied a larger counts 
limit; these will be mentioned explicitly. In some instances 
entire observations did not contain enough counts to reach the 
16,000 counts minimum; in those cases we created a single 
data point from the whole observation if the total number of 
counts was larger than 10,000, but excluded it from further 
consideration if it had fewer counts. (Any exceptions to this 
will be explicitly mentioned.) We note that for most observa¬ 
tions of the four bright sources studied here no rebinning was 
necessary, since each 16 s time bin usually contained more 
than the minimum 16,000 counts. With this minimum counts 
limit the uncertainty due to counting statistics is at most ^2% 
in the soft color and at most ^3% in the hard color (often 
much less); the uncertainty in the intensity is always less than 
1 %. 


3. RESULTS 
3.1. XTE J1701-462 


XTEJ1701-462 was discovered with the RXTE All-Sky 
Monitor (ASM) on 2006 January 18 as the source was en¬ 
tering outburst ( [Remillard & Lin|2006| . The first RXTE PCA 
observation took place on 2006 January 19, and the source 
was subsequently observed with the PCA on average ~ 1.5—2 
times per day for the remainder of the 19 month outburst, 
apart from a ~50 day period in late 2006 and early 2007 dur¬ 
ing which the source could not be observed due to proximity 
to the Sun. The source returned to quiescence in early 2007 
August and has since remained inactive apart from occasional 
low-level flaring (up to ~10 35 erg s -1 ), which has been ob¬ 
served by S wift an d XMM-Newton (Fridriksson et al. 2010 


201 1| ). As in |H10| dates during the outburst will in this paper 


be referred to as days since the start of 2006 January 19 (MJD 
53754.0). 


The combined CD and HID for the entire 2006-2007 out¬ 
burst of XTE J1701-462 are shown in Figure [l] the diagrams 
clearly illustrate that the source tracks exhibited strong sec¬ 
ular motion during the outburst. As in |H10| we divided the 
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Figure 2. CDs and HIDs for XTE J1701-462, showing the secular evolution of the source. Each of the 12 panels (consisting of a CD on the left and an HID 
on the right) corresponds to a particular selection of data from the entire RXTE PC A data set for the source; red data points correspond to a particular subset 
(exhibiting minimal secular shifts) within a given selection (see Table [2}. Particular source states/branches are indicated in a few of the panels (see text for 
definitions). The dashed and solid lines show the approximate paths followed by the NB/FB and HB/NB vertices, respectively, as the tracks evolve. Data are 
binned to a minimum of 16,000 counts per data point, except for the IS and EIS data in selection L, which are binned to a minimum of 64,000 counts. Note 
the change in the intensity scale between the different HID columns. We also note that essentially the same tracks (with only small modifications here) were 
presented in slightly different form in Figures 3 and 4 in |H10| 


observations of the source into 12 subsets—hereafter referred 
to as selections—which we label A-L. In Figure |2] we show 
the CDs and HIDs for all 12 selections, which illustrate the 
secular evolution of the source during the outburst. We la¬ 
bel the source states/branches of a few representative tracks: 
a Cyg-like Z track (selection B), a close-to Sco-like Z track 
(E), and two atoll-like tracks (K and L). For the Z-like tracks 
these are the horizontal branch (HB), the HB upturn, the nor¬ 
mal branch (NB), the flaring branch (FB), and the dipping FB. 
For the atoll-like tracks these are the lower and upper banana 
branches (LB and UB), the island state (IS), and the extreme 
island state (EIS). The banana branch, IS, and EIS are also re¬ 
ferred to as the soft, intermediate/transitional, and hard state, 


respectively. (We refer to H10 for examples of CDs/HIDs of 
several Cyg/Sco-like Z sources and atoll sources and a com¬ 
parison of those with the CD/HID tracks of XTE J1701-462.) 

The 12 data selections mostly correspond to particular 
ranges in the low-energy (2.0-2.9 keV) count rate, which [HT0| 
found to closely trace the secular changes during most of the 
outburst. This tracing of the secular evolution with the low- 
energy count rate breaks down when the source enters the HB 
or a dipping FB. For days 0-28 time-based selections (A and 
B) were therefore used, rather than ones based on the low- 
energy count rate. In addition, data from several HB excur¬ 
sions taking place after day 28 had to be moved to the same se¬ 


lection as neighboring NB data in the low-energy light curve. 
The time and count rate intervals we used are given in Ta¬ 
ble J2 ] Figure 2 in |H10 shows a low-energy light curve of the 
entire outburst, indicating the count rate intervals used there 
(which are for the most part equivalent to ours) and the HB 
data moved between selections. 

Within most of the selections some secular motion is still 
evident, mostly in the HID. We therefore in 10 of the 12 se¬ 
lections color red the data points from a subset of the ob¬ 
servations used in each case, to show tracks with minimal 
secular shifts; the subintervals in time or count rate defining 
these subsets are given in Table [2] We also show, in both 
the CDs and HIDs, two straight lines that the NB/FB (lower) 
and HB/NB (upper) vertices of the tracks are seen to follow 
closely (as pointed out by Lin et al. |2009b| . We note that se¬ 
lections A-L do not represent a strict monotonic progression 
in time (although overall the evolution in time was from A to 
L); the source moved back and forth between selections dur¬ 
ing the outburst, mostly within the range F-H (see Figure 2 
in |H10| ). It is clear from the sequence in Figure [2] that going 
from panel A to L the tracks smoothly evolve in shape from 
Cyg-like Z to Sco-like Z to atoll tracks as the upper and lower 
vertices both move monotonically along the vertex lines to 
higher color values and the overall intensity of the tracks de¬ 
creases. 
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Table 2 

Time or Count Rate Intervals Used for Data Selections 


Selection 

Full Interval 
(counts s“‘ PCLT 1 ) 3 

Subinterval 
(counts PCIT 1 ) 3 

A 

Days 13.5-20.5 

Days 13.5-19.5 

B 

Days 0-13.5 and 20.5-28 

Days 2-13.5 

C 

106-114/Days 28-32.5 


D 

83-106 

89-97 

E 

66.5-83 

72.5-78 

F 

56-66.5 

56-58.36 

G 

46-56 

50-52 

H 

38.5-46 

39.6-43 

I 

31-38.5 

33.5-35.5 

J 

23.5-31 

26-28.5 

K 

15-23.5 

20-22.5 

L 

0.4-15/Days 550-564 



Notes. 

a Count rates are in the 2.0-2.9 keV band. Time intervals refer to days since 
MJD 53754.0. 


After day 566 (as the outburst was ending) the 
count rate from the source reached a roughly con¬ 
stant level of ^2 counts s -1 PCU -1 in the 2-60 keV band 
(^0.2 counts s -1 PCU -1 in the low-energy band). This resid¬ 
ual flux can be attributed to diffuse Galactic background emis¬ 
sion (see |Fridriksson et al.|2010| . Figure [3] shows an alterna¬ 
tive version of panel L in Figure [2] where we have binned the 
data more heavily and have subtracted this diffuse background 
emission based on data from observations made during the 
three weeks following the end of the outburst. In addition, we 
include in this version data from two observations made dur¬ 
ing days 564-566 (constituting the leftmost data point), which 
had a low-energy count rate of ~0.36 counts s -1 PCU -1 . We 
furthermore use a logarithmic scale for the horizontal axis in 
the HID so as to better illustrate the behavior of the source at 
the lowest count rates. 


3.2. CygX-2 

Cygnus X-2 (Cyg X-2) is one of the longest-known and 
most extensively studied X-ray binaries. Cyg X-2 was clas¬ 
sified as a Z source by Hasinger & van der Klis ( |1989j ) based 
on EXOSAT data and is the prototype of the Cyg-like sub¬ 
group of the persistent Z sources. However, as was clear al¬ 
ready from pr e-RXTE data Cyg X-2 is unique among the six 
“classic” (persistent Galactic) Z sources in that it shows by 
far the strongest secular evolution (e.g., |Kuulkers et al.|1996 
Wijnands et al. ]T997j>. 


3.2.1. Analysis 

In Figure [4] we show a CD/HID based on the entire RXTE 
PC A data set for the source; strong variations in the shape 
and location of tracks are readily apparent. This combined 
CD/HID for Cyg X-2 has strong similarities to the one for 
XTEJ1701-462 (Figure [TJ). However, the fragmented na¬ 
ture of the data set (obtained over a period of 15 years) 
forces us to analyze the data in a manner different from the 
XTE J1701-462 analysis. 

For this analysis we considered all 591 individual ObsIDs 
for pointed observations of Cyg X-2 made during the RXTE 
mission (13 of which contained no useful data). Often mul¬ 
tiple observations were fairly densely clustered together in 



Figure 3. Alternative version of panel L in Figure|2jwith the data binned into 
fewer groups and Galactic background emission subtracted. Data in the LB 
(i.e, intensity >65 counts s -1 PCU -1 and hard color <0.85) are binned to a 
minimum of 128,000 counts per data point; data in the IS and EIS are binned 
to approximately one data point per day. The data represented by the leftmost 
point were not included in Figure [2] 

time. Going chronologically through the data set the source 
jumps erratically back and forth around the CD/HID, in any 
given observation usually tracing out only short partial track 
segments. As a first step we went through the entire data set 
in time order, at any point combining into a single CD/HID 
track as many consecutive ObsIDs as possible without intro¬ 
ducing clearly noticable secular shifts. We refer to such clus¬ 
ters of data (containing all data from a certain time interval) as 
subsets. In a few cases significant secular motion took place 
during a single ObsID, requiring the observation to be split 
up between two different subsets. To identify secular motion 
we examined the tracks not only in the CD and HID, but also 
in a diagram of soft color versus intensity (which we refer to 
as an SID). The SID frequently yielded useful extra informa¬ 
tion/constraints. The result of this process was ^300 subsets, 
which span periods as long as several days, although the vast 
majority is less than a day in length. 

Considering only individual subsets results in a few tracks 
that seem to a large extent complete, but most of the subsets 
consist of shorter segments (i.e., incomplete tracks). By an 
“incomplete” track we mean that had the source stayed long 
enough at a particular stage in its secular progression (and 
had our observational coverage of the source been sufficiently 
comprehensive during that time) we expect that the source 
would have traced out a fuller (i.e., continuous and possibly 
more extended) track. It is important to note that this “com¬ 
pleteness” is a function of location in the secular progression, 
as not all branches of the Z track are present at every stage 
of the secular evolution, and that there is always some uncer¬ 
tainty regarding what constitutes a “complete” track. 



1.2 1.4 1.6 1.8 2 2.2 2.4 500 1000 1500 2000 2500 

Soft color Intensity (counts s ” 1 PCU” 1 ) 

Figure 4. CD/HID representing the entire RXTE PCA data set of Cyg X-2. 
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Figure 5. A sequence of CDs and HIDs for Cyg X-2 illustrating the secular evolution of the source. Table [3] lists the data used in each panel. In general, a 
panel combines data from observations widely separated in time (by months or years); exceptions to this are panels A, B, K, and N. In each of the other panels, 
a particular subset of the data—obtained over a period of at most a few days—is shown in red. The dashed and solid lines show the approximate paths followed 
by the lower and upper vertices, respectively. Note the change in the intensity scale between the different HID columns. 


Fortunately, it was in general possible to combine subsets 
from various times throughout the RXTE mission to form 
more complete tracks than are available from individual sub¬ 
sets. Usually, the most complete individual-subset tracks 
served as a foundation on which these combined tracks were 
built. To guide us in this process we took advantage of over¬ 
lapping track segments and made sure that they lined up in 
all three diagrams—the CD, HID, and SID—which together 
provided a fairly stringent criterion for the appropriateness of 
combining particular segments. 

In Figure [5] we show a sequence of 14 tracks, chosen (from 
a larger set of tracks) to illustrate as clearly as possible the 
overall secular evolution exhibited by the source—both of the 
individual branches and the locations of the tracks—as they 
move through the CD/HID. As is the case for XTE J1701-462 
the lower and upper vertices approximately follow straight 
line paths in both the CD and HID as the tracks evolve and 
shift in the diagrams. We show illustrative lines in Figure [5] 
and order the panels based on the vertex locations of the 


tracks, starting (panel A) at the highest intensities and low¬ 
est color values. 

The number of individual subsets in the tracks in the 
Cyg X-2 sequence ranges from 1 to 10. The subsets (i.e., 
time intervals) and ObsIDs used for each panel in Figure [5] 
are listed in Table [3] For the panels that consist of more than 
one data subset (i.e., all except A, B, K, and N), a single repre¬ 
sentative subset is plotted in red. The combined exposure time 
of the data shown in Figure [5] is ~31% of the total exposure 
time of the Cyg X-2 data set (see Table [l}. As far as we can 
tell practically all the remaining data seem to belong to tracks 
similar to (and usually intermediate between) the ones shown 
in Figure [5] (i.e., none of the other subsets is inconsistent with 
belonging to such tracks). 

Although the tracks in panels K and L do not line up well 
enough in the HID (contrary to the CD) to be combined into a 
single track, together they should give a reasonably good idea 
of the rough shape of a (near-)complete track at this point in 
the secular progression. It is clear that such a track is similar 
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Table 3 

Time Intervals and Observations Used to Create Cyg X-2 Tracks 




Interval Length 


Exp. Time 

Panel 

Time Interval (MJD) 

days (hr) a 

ObsIDs 

ks (hr) b 

A 

50316.592-50316.764 

0.17 (4.1) 

10063-[09:10J-01-00 

9.0 (2.5) 

B 

51697.550-51698.705 

1.15 (27.7) 

40019-04-05- [00:10] 

33.9 (9.4) 

C 

50996.450-51000.593 (R) 

4.14(99.4) 

30418-01-[(01:05)-00,02-01] 

37.4 (10.4) 


51561.344-51561.517 

0.17 (4.1) 

40017-02-19-00 

10.1 (2.8) 


53138.767-53138.855 

0.09 (2.1) 

90030-01-16-00 

5.2 (1.4) 


53788.496-53788.577 

0.08 (2.0) 

91009-01 -42- [00:01 ] 

2.7 (0.7) 


54649.890-54650.626 

0.74 (17.6) 

93443-01-01-[02,14:21] 

23.6 (6.6) 

D 

51009.488-51009.725 

0.24 (5.7) 

30046-01-01-00 

12.9 (3.6) 


51349.631-51349.669 

0.04 (0.9) 

40017-02-09-01 

3.3 (0.9) 


53079.687-53079.785 

0.10(2.3) 

90030-01-04-00 

5.2 (1.4) 


53222.440-53222.532 

0.09 (2.2) 

90030-01-33-[00:01] 

4.8 (1.3) 


53286.420-53286.512 

0.09 (2.2) 

90030-01-46-00 

5.3 (1.5) 


53291.340-53291.433 

0.09 (2.2) 

90030-01-47-00 

5.4 (1.5) 


54007.232-54007.331 

0.10(2.4) 

92039-01-15-00 

5.7 (1.6) 


54009.195-54010.341 

1.15 (27.5) 

92039-01-[17,18]-00 

8.2 (2.3) 


54452.324-54452.554 

0.23 (5.5) 

90022-08-04- [00:01 ] 

9.5 (2.6) 


54648.066-54649.329 (R) 

1.26 (30.3) 

93443 . 0 1-01- [00,000,03:13] 

37.2 (10.3) 

E 

50168.900-50169.731 (R) 

0.83 (20.0) 

10066-01-01- [00,000,001] 

37.3 (10.4) 


51048.360-51048.599 

0.24 (5.7) 

30046-01-07-00 

13.4 (3.7) 

F 

51055.495-51055.782 

0.29 (6.9) 

30046-01-08-00 

13.4 (3.7) 


51266.719-51267.126 

0.41 (9.8) 

40017-02-05- [00:01 ] 

15.9 (4.4) 


51413.523-51413.746 

0.22 (5.4) 

40017-02-12-00 

13.0 (3.6) 


52429.372-52429.908 

0.54 (12.9) 

70016-01-01-[02,04:06] 

12.4 (3.4) 


53511.868-53512.546 

0.68 (16.3) 

91010-01-01-[01,06:09] 

29.9 (8.3) 


53685.118-53685.193 

0.08 (1.8) 

90030-01-82-00 

4.3 (1.2) 


54395.861-54398.037 (R) 

2.18(52.2) 

92038-01-[08-(00:02),09-(00:02),10-(00:03)] 

12.3 (3.4) 

G 

51081.447-51081.725 

0.28 (6.7) 

30046-01-12-00 

14.3 (4.0) 


51444.401-51445.215 

0.81 (19.5) 

40019-04- [01 -01 G,01 -02G,01 -03,02-00,02-000] 

17.6 (4.9) 


52533.929-52535.087 

1.16(27.8) 

70015-02-[01-01G,02-00] 

19.1 (5.3) 


53311.026-53311.105 

0.08 (1.9) 

90030-01-51-[00:01] 

2.6 (0.7) 


53335.431-53335.724 

0.29 (7.0) 

90022-08-01 -00, 90030-01 -56-[00:01 ] 

8.2 (2.3) 


53375.032-53375.053 

0.02 (0.5) 

90030-01-64-00 

1.8 (0.5) 


53394.710-53394.812 

0.10(2.5) 

90030-01-68-00 

6.3 (1.8) 


53488.223-53488.310 

0.09 (2.1) 

91009-01-12-[00:01] 

3.4 (0.9) 


53493.152-53493.227 (R) 

0.07 (1.8) 

91009-01 -13- [00:01 ] 

2.8 (0.8) 


53498.069-53498.142 

0.07(1.8) 

91009-01-14-00 

3.1 (0.9) 

H 

51539.192-51539.562 (R) 

0.37 (8.9) 

40017-02-18-[00:02] 

14.4 (4.0) 


53119.069-53119.174 

0.10(2.5) 

90030-01-12-00 

6.8 (1.9) 


53330.721-53330.805 

0.08 (2.0) 

90030-01-55-[00:01] 

2.9 (0.8) 


53347.497-53347.527 

0.03 (0.7) 

90022-08-02-00 

2.6 (0.7) 


53744.187-53744.278 

0.09 (2.2) 

90030-01-94-00 

5.2 (1.4) 


54388.862-54389.040 

0.18(4.3) 

9203 8-01-01- [00:02] 

3.8 (1.1) 

I 

50629.832-50631.933 (R) 

2.10(50.4) 

20053-04-01-[00:04,06,010,020,030] c 

77.5 (21.5) 


51029.751-51030.032 

0.28 (6.7) 

30046-01 -04- [00:01 ] 

13.8 (3.8) 


53266.795-53266.821 

0.03 (0.6) 

90030-01-42-00 

2.3 (0.6) 


53758.990-53759.021 

0.03 (0.7) 

90030-01-97-00 

2.7 (0.7) 


54427.903-54428.970 

1.07 (25.6) 

93082-02-04-[00:01] 

2.4 (0.7) 

J 

53123.996-53129.015 (R) 

5.02(120.5) 

90030-01 - [ 13:14] -00 

8.8 (2.4) 


53320.893-53320.934 

0.04(1.0) 

90030-01-53-[00:01] 

1.4 (0.4) 

K 

51528.887-51529.451 

0.56(13.5) 

40019-04- [03-01,04-(00:01),04-000] 

22.8 (6.3) 

L 

50719.393-50719.759 

0.37 (8.8) 

20057-01-01- [00,000] 

20.2 (5.6) 


51536.473-51536.968 (R) 

0.50(11.9) 

40021-01-02- [00:04] d 

19.2 (5.3) 

M 

51535.886-51536.436 (R) 

0.55 (13.2) 

40021-01- [01-01,02-00,02-000] d 

26.7 (7.4) 


53202.754-53202.834 

0.08 (1.9) 

90030-01-29-00 

4.3 (1.2) 

N 

51535.010-51535.571 

0.56(13.5) 

40021-01-01- [00,000,02] 

25.8 (7.2) 


Notes. Subsets colored in red in Figure[5]are denoted by (R) in the Time Interval column. In the ObsIDs column a colon denotes a range. 
a The interval length is shown in units of both days and hours. 
b The total exposure time is shown in units of both ks and hours. 

c As the source showed significant secular motion during observation 20053-04-01-04, only an early part of it was used. 

d The source showed clear secular motion during observation 40021-01-02-00. Data from the first two orbits were used in panel M and data from the third (and final) orbit was used in 
panel L. 
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in shape to those of the Sco-like Z sources. 

We note that small secular shifts do occur in some of these 
tracks; some of the individual segments used have a broad ap¬ 
pearance (e.g., the red-colored segment in panel D), which is 
in some cases probably due to mild secular motion, and in 
some instances we matched up segments despite their not lin¬ 
ing up perfectly in the HID or SID if the overall appearance 
of the track is only minimally affected by this. There was 
inevitably sometimes some ambiguity regarding whether cer¬ 
tain segments used in a given track rather belonged to a track 
slightly shifted from the one in question. In particular, the 
HB upturn was usually observed in short isolated segments 
in the CD/HID, and it was often hard to judge exactly with 
which tracks those should be combined. However, the overall 
conclusions about the secular evolution of the source that we 
infer from the data are not sensitive to these ambiguities in the 
combining process. Finally, we note that many of the tracks in 
Figure [5] likely still suffer from some incompleteness, either 
because data segments that would serve to complete them are 
simply not available in the Cyg X-2 data set, or because we 
felt that there was too much ambiguity in whether candidate 
segments were appropriate. In particular, we conclude that the 
track in panel A is likely missing the HB, the tracks in panels 
B and F are presumably missing (most of) the HB upturn, the 
one in panel H is missing most of the HB, the track in panel J 
is missing the NB and HB, and the one in panel K is missing 
theFB. 


3.2.2. Comparison with XTE J1701-462 

Cyg X-2 exhibits secular evolution that is for the most part 
very similar to that of XTE J1701-462. As the overall in¬ 
tensity decreases the tracks smoothly evolve in shape from 
Cyg-like to Sco-like Z tracks. The tracks in Cyg X-2 panels 
A-E (Figure [5]) are very similar to the Cyg-like Z tracks of 
XTE J1701-462 in panels A and B. Taken together, the tracks 
in Cyg X-2 panels K and L show that at this point in the sec¬ 
ular progression of the source its tracks look similar to (and 
perhaps intermediate between) those in XTE J1701-462 pan¬ 
els E and F, the latter of which is very similar to the tracks of 
the persistent Sco-like Z sources, such as GX 17+2 (see, e.g., 
H10l|Lin et al.|2012|). Finally, the tracks in Cyg X-2 panels M 
and N look very similar to the tracks in XTE J1701-462 pan¬ 
els G and H. As the Cyg X-2 tracks shift and evolve in shape, 
both the upper and lower vertices quite closely follow straight 
lines in the CD and HID; as the overall intensity of the tracks 
decreases, the NB is squeezed between the converging vertex 
lines and gradually shortens, as observed in XTE J1701-462. 

The main differences between Cyg X-2 and 
XTEJ1701-462 lie in the fact that Cyg X-2 is persis¬ 
tently luminous and has never been seen to enter the atoll 
regime (although the track in panel N does not seem to have 
an NB, short NB segments seem to exist down to the very 
lowest intensities observed in the Cyg X-2 data set) and in 
the somewhat different FB behavior exhibited by Cyg X-2. 
In both sources the FB gradually evolves from being a purely 
dipping FB in the earliest (highest-intensity, Cyg-like) panels 
to a “proper” flaring (Sco-like) FB in later panels. However, 
in Cyg X-2 the rotation of the FB is for the most part in the 
opposite direction to that seen in XTE J1701-462 in both the 
CD and HID. In Cyg X-2 the FB also develops a more com¬ 
plicated morphology in the HID, making a counterclockwise 
twist of between 180° and 270° around the lower vertex in 
panels G and H (such a twistAvas also observed in EXOSAT 
data by Kuul kers et al. 1996| ) and then assuming a jagged 


S-like shape in panels I and J that gradually straightens out 
going from panel I to N. We note that a similar, although less 
pronounced, S-shaped FB was observed in the HID of panel 
E in XTE J1701-462. 


3.3. CirX-1 

Circinus X-l (CirX-1) features some of the richest and 
most complex phenomenology seen among the known neu¬ 
tron star X-ray binaries, and although it has been exten¬ 
sively studied for over four decades many of its properties 
remain poorly understood. The binary has an orbital period 
of ~16.6 days, ide ntified from periodic fla ring first observed 
in the X-ray band (Kaluzienski et al.||1976|), and later in the 
radio, infrared, and optical ( |Whelan et al. 11977} |Glass|197~8| 
Mon eti|1992| ). These flares are thought to be due to enhanced 


accre tion near periastron passage in a highly eccentric orbit 


1 P ; 

(e.g., Murdin et al. 1980; Tauris et al 


2004). Heinz et al. (20 13} identify the radio nebula sur- 
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rounding the binary as a relatively young (<5000 yr old) su 
pernova remnant, which makes Cir X-l the youngest known 
X-ray binary and provides an explanation for the eccentricity 
of the orbit. This implies that the donor is likely an early-type 
star (see also Jonker et al. 2007 ). However, type I X-ray bursts 
detect ed from the source ( Tennant et al.|1986bla||Linares et al.| 
2010| ) show that the neutron star is weakly magnetized, which 
is very unusual for a high-mass companion. 

The X-ray emission from CirX-1 is highly variable on a 
wide range of timescales. In Figure[6]we show an RXTE ASM 
light curve of the source. The source was in a historic ally high 
state during the first few years of the RXTE mission ( [Parkin¬ 
son et al. |2003| ), with an average flux of ~1.3 Crab (and a 
maximum observed flux with the ASM of ^3.5-4 Crab). The 
flux started gradually decreasing in mid-to-late 1999, and kept 
doing so until the source became undetectable with the ASM 
and showed no measurable activity over a two-year period in 
2008-2010. In 2009 the source was observed with Chandra 
at a flux of only a few tenths of a milliCrab ( [Sell et al.|2010| ). 
Since 2010 May, however, the source has shownsporadic ac¬ 
tivity. When active the source usually exhibits complex vari¬ 
ability over the course of an orbital period, featuring both ab¬ 
sorption dips and flaring in the X-ray band; the effects of this 
on the ASM light curve can easily be seen in Figure [6] In 
the inset we show five orbital cycles from the bright phase 
of the source early in the RXTE mission. X-ray flaring is 
strongest during the first few days after it commences; it then 
gradually decreases in strength as the orbit progresses. The 
source usually exhibits strong X-ray dipping during the last 
^0.5—1 days before the onset of flaring and then intermittent 
dips for up to two days afterwards ( [Shirey et al.||1998| ); we 
di scuss this dipping in more detail below! 

|Oosterbroek et al.| (|1995|) reported indications of both Z 
and atoll behavior in EXOSAT data of CirX-1. |Shirey et al. 
( |1999a| ) analyzed data from 10 days of RXTE PCA ob- 
servations in 1997. They detected all three Z branches— 
although with some differences in shape compared to the clas¬ 
sic Z sources—and found that the tracks moved around the 
CD/HID and evolved in shape. The identifications of the Z 
branches were supported by timing data. 


3.3.1. Analysis 

In our analysis we considered all 811 individual ObsIDs for 
pointed observations of Cir X-1 made during the RXTE mis¬ 
sion (18 of which contained no useful data). A CD/HID using 
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Figure 6. RXTE ASM light curve of Cir X-1 covering the entire lifetime 
of the RXTE mission. Data points in the main plot are one day averages. 
The inset zooms in on five orbital cycles in early 1997; each data point there 
corresponds to a single dwell of the ASM. The long-dashed lines in the main 
plot and inset show the typical ASM count rate level for the Crab Nebula. 


Time (year) 


2000 2005 2010 



the entire data set (which spans ~15 years) is shown in the 
upper plot in Figure [ 7 ] The diagrams are heavily affected by 
both absorption and secular shifts and shape changes in the 
source tracks. As mentioned above, most of the dipping oc¬ 
curs close to the time of presumed periastron pa ssage, which 
is often associated with flaring. |Clarkson et al.| ( 2004] ) fitted 
a quadratic ephemeris to the times of dips observed in RXTE 
ASM data from 1996 to 2003 and found that it provides a 
good predictor of the X-ray light curve. The strongest dip¬ 
ping often produces a characteristic track in the CD and HID; 
this track goes up to high color values and has two sharp bends 
in the CD. This can be seen in Figures [ 7 ] and [8} in the latter 
we show various diagrams for 7 days of observations in 1997 
June, which we discuss further below. (The observations used 



-ray emission is composed of two components 
a bright component that is subject to heavy absorption and a 
faint component that is unaffected by the absorption, the lat¬ 
ter perhaps due to X-rays from the central source scattered 
into our line of sight by surrounding material (Shirey et al. 


1 1999b) ). We note that the source also often exhibits shallower 
dips (also seen in Figure[8]) which do not result in a (full) track 
of the sort described above. 

Before proceeding with our analysis of the Cir X-l data set 
we removed by eye (to the extent possible) data points af¬ 
fected by absorption dips. An example of this is shown in 
Figure [8j where removed data points are colored red. As can 
be seen in the figure, during these observations the source 
showed intense dipping shortly before the start of flaring (at 
day ^611.5), producing the characteristic dipping track in the 
CD/HID. Some shallower dips were then observed during the 
first days of flaring. We note that the periodic “flaring” 
during an orbital cycle is in general associated with motion 
along all three (Cyg-like) Z branches, rather than being exclu¬ 
sively associated with motion along the FB (which is mostly 
a dipping FB in these particular observations). As is apparent 
from the figure, it is almost impossible to identify absorption 
dips on the basis of the light curve alone during periods of 
flaring. However, tracking the behavior of the soft color (and, 
to a lesser extent, the hard color) as a function of time, as 
well as inspecting the CD and HID, greatly aids in identify¬ 
ing dipping. In addition, we took into account the dipping 
ephemeris of|Clarkson et al.|(|2004|) when performing the re¬ 
moval, since the vast majority of dipping events take place in 
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Figure 7. CD and HID representing the entire RXTE PC A data set of Cir X-l 
before (upper plot) and after (lower plot) the removal of data affected by 
absorption. A small portion of the data at very high soft and hard color values 
falls outside the diagrams in the upper panel. Note the differences in scale in 
soft and hard color between the upper and lower panels. 

the ~ 1—1.5 days immediately before or after phase 0. How¬ 
ever, we note that (shallow) dips are sometimes seen later in 
an orbital cycle; one example can be seen in Figure [8] In 
general, when unsure whether a given data segment was af¬ 
flicted by dipping or not, we in general opted to rather err 
on the side of caution and remove the segment in question. 
The “cleaned” CD/HID resulting from our manual removal 
of absorption-affected data is shown in the lower plot in Fig¬ 
ure |71 It was of course unavoidable that some instances of 
minor dipping remain and that a small amount of unaffected 
data be removed. However, we expect that the effects of this 
on the conclusions we draw from the data are negligible. 

After removing data points affected by absorption dips, as 
described above, we organized the data into subsets, similar to 
our Cyg X-2 analysis described in the previous section. This 
resulted in ^300 subsets spanning periods as long as several 
days, although most are less than a day in length. From these 
we created a sequence of 12 CD/HID (partial) tracks, shown 
in Figure[9] These tracks were chosen to illustrate as best pos¬ 
sible the overall secular behavior of the source. In all cases 
except two the data in a given panel are from a single sub¬ 
set, whose time intervals range from ^2 hr to ~7 days. The 
exceptions are panels A and I, where 2-3 independent seg¬ 
ments (shown in different colors) were combined. It was not 
viable to form a whole sequence of combined and more com¬ 
plete tracks, as we did for Cyg X-2, since the data for Cir X-l 
show a much greater dynamic range and do not sample the 
HID densely enough for this. Overall, the upper and lower 
vertices of the Cir X-l tracks move systematically up and to 
the right in the CD, and up and to the left in the HID, as 
the tracks evolve in shape—similar to the behavior seen in 
XTE J1701-462 and Cyg X-2. The paths of the vertices can 
be approximated by the straight lines shown in Figure [9] We 
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Figure 8. Diagrams demonstrating the removal of data affected by absorption from 7 days of observations of Cir X-l in 1997 June; the left plot shows intensity 
and soft color as a function of time, and the right plot shows a CD/HID. We identified the red data points with absorption and removed them from the dat a set 
before proceeding with further analysis. In the left plot the dashed vertical line shows the time of zero phase according to the dipping ephemeris oflClarkson et al.l 
{2004}. 


order the tracks based on the position of the lower and/or up¬ 
per vertex in the CD. (When both are present, the two vertices 
give consistent results.) In cases where some ambiguity re¬ 
mains due to very similar vertex locations (in particular, pan¬ 
els F and G), we use the shape of the track in the CD to decide 
the ordering, in particular the FB (i.e., we choose the ordering 
that produces a gradual evolution in the shape of the track in 
the CD). The combined exposure time of the data shown in 
Figure [9] is ~10% of the total exposure time of the Cir X-l 
data set. The data subsets (time intervals) and ObsIDs used in 
each panel are listed in Table [4] 

We note that a few of the tracks in Figure[9]show some signs 
of secular motion in the HID (but much less in the CD); how¬ 
ever, in most cases this has a negligible effect on the overall 
appearance of the track. The instance where this is most not- 
icable is in the HID of panel E, where the upper parts of the 
track (NB and especially HB) show shifting toward lower in¬ 
tensities relative to the FB. A similar, but smaller, shift affects 
the HID in panel F. We also note that the behavior of Cir X-1 
is not as regular as that of XTE J1701-462 and Cyg X-2, and 
the vertices in the HID (and to a lesser extent in the CD) often 
deviate significantly from the lines shown. As will be dis¬ 
cussed in Section [4] the relationship between a track’s shape 
and its location in the HID is also not nearly as tight for 
Cir X-l as for XTE J1701-462 and Cyg X-2, where the mo¬ 
tion of both vertices in both the CD and HID, as well as the 
decrease in overall intensity, is monotonic along the entire se¬ 
quence. The behavior of Cir X-l is more regular in the CD 
than the HID, which is why we chose to use the CD as the ba¬ 
sis for our ordering of the tracks. However, apart from some 
irregularities in panels F-I, the Cir X-l sequence does show a 


gradual decrease in overall intensity. 

Panel L shows what seems to be an atoll transition from 
the banana branch through the IS and to an EIS. These data 
were obtained in 2010 May/June as the intensity steadily de¬ 
creased over a ~7 day period (part of a ^50 day minioutburst 
following two years of nondetection by the RXTE ASM). Ob¬ 
servations preceding this 7 day period show what looks like a 
sequence of segments at successively lower count rates from 
tracks similar to the one in panel K. Figure [lO] shows the same 
data as in panel L (in black) but with more binning. D’Ai et al. 
( |2012| ) studied the spectral evolution of the source during the 
entire 2010 May-June minioutburst with spectral fitting and 
also concluded that the source transitioned from the atoll soft 
to hard state. In 2010 August the source showed a similar tran¬ 
sition (again with an intensity decline over ^7 days), which 
we show with red data points in Figure [lO] We note that the 
apparent decrease in hard color with decreasing intensity in 
the EIS in both cases may well be due to soft diffuse back¬ 
ground emission (which could not be subtracted) affecting the 
data points at the lowest intensities. 

3.3.2. Comparison with XTE J1701^62 and Cyg X-2 

Like XTE J1701-462, Cir X-l has been observed in all NS- 
LMXB subclasses (Cyg-like Z, Sco-like Z, atoll). The over¬ 
all secular evolution has many similarities to that observed 
for XTE J1701-462 and Cyg X-2, and many of the individual 
tracks have shapes similar to those seen for those two sources. 
Progressing along the Cir X-l sequence the NB grows shorter 
and rotates counterclockwise in the HID, while the HB ro¬ 
tates clockwise as it shortens, similar to what was seen for 
XTE J1701-462 and Cyg X-2. As in those two sources the 
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Figure 9. A sequence of CDs/HIDs for Cir X-l illustrating the secular evolution of the source. The data points in each panel were obtained within a relatively 
short period of time, ranging from ~4 hr to ~7 days, except for panels A and I, which combine three and two segments, respectively, widely separated in time 
(with each segment shown in a particular color). The dashed and solid lines are lower and upper vertex lines similar to those shown for XTE J1701-462 and 
Cyg X-2. Data are binned to a minimum of 16,000 counts per data point, except for IS and EIS points in panel L, which are binned to a minimum of 32,000 counts. 
Note the change in the intensity scale between the different HID columns. 


FB shows the most complex behavior of the three branches 
in Cir X-l. In the CD of Cir X-l the FB gradually rotates 
counterclockwise and evolves in shape in a similar fashion to 
XTE J1701-462. In the HID the behavior of the FB is more 
irregular in Cir X-l than XTE J1701-462 and Cyg X-2, but 
overall it seems to rotate counterclockwise as it evolves from 
a dipping FB in the higher-intensity tracks to a Sco-like FB at 
lower intensities. 

In panel G, the Cir X-l FB in both the CD and HID has 
a shape similar to those of the pe rsistent Cyg-like Z sources 
GX 340+0 and G X 5-1 (see, e.g., |Jonker et al.|1998[ |Jonker 
|et al.|2002l|H10| ). We note that in the CD of XTE J 1701^62 


the shape and orientation of the FB in panels A and B seems 
intermediate between Cir X-l panels E and F, whereas the 
XTE J1701-462 FB in panels C (which has a very incomplete 
track) and D seems intermediate between Cir X-1 panels G 
and H. XTE J1701-462 may therefore have traced out a CD 
(and HID) track similar to those of GX 340+0 and GX 5-1 in 
between panels B and C; this portion of the secular progres¬ 
sion of the source was missed due to a gap in RXTE coverage 
and rapid secular evolution. However, we note that in none of 
the four sources compared in this paper do we see a full track 
in both the CD and HID where the shapes of all three spectral 
branches closely match those of GX 5-1 and GX 340+0. 

The atoll transitions of Cir X-l—especially the one in 2010 
May/June—resemble the one of XTE J1701-462 (see Fig¬ 


ures [3] and [T0|; one notable similarity is that the soft color 
decreases significantly throughout the transitions, in contrast 


to the increase observed in most atoll sources (e.g., [Fridnk s- 
2011). However, there is also a striking difference be- 


son 


tween the atoll transitions of the two sources: the data from 
the banana branch preceding the ascent to the hard state ex¬ 
tend to very high soft color for Cir X-l, occupying parame¬ 
ter space in the CD never explored before by the source in 
the RXTE archive except during absorption dips or (partly) 
during traversals to the tip of the (nondipping) FB. The data 
points in question are unlikely to be associated with the FB, 
given the small changes in intensity and hard color observed, 
and the long timescale involved (more than a day). This 
pre-atoll-transition behavior of Cir X-l is in stark contrast to 
XTE J1701-462, where no such excursion to high soft color 
values was seen, and the region occupied by the source in 
the CD in selection L before moving to the IS was a logical 
extension of the movement of the source in the preceding se¬ 
lections. We also note that the observations in 2010 May/June 
were the first time in almost 15 years of RXTE observations 
that Cir X-l was seen to transition to the atoll hard state. 

The Cir X-1 track in panel D bears a strong resemblance 
in shape to the highest-intensity tracks (panels A and B) of 
XTE J1701-462 (in particular when the HID is plotted on the 
same scale) and occupies a similar color range (especially in 
hard color). But Cir X-l also shows tracks that reach much 
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Table 4 

Time Intervals and Observations Used to Create Cir X-l Tracks 




Interval Length 


Exp. Time 

Panel 

Time Interval (MJD) 

days (hr) a 

ObsIDs 

ks (hr) b 

A 

50365.216-50365.307 

0.09 (2.2) 

10068-08-02-00 

5.5 (1.5) 


50497.364-50497.437 (G) 

0.07(1.8) 

20095-01-01-00 

4.8 (1.3) 


50711.692-50711.763 (R) 

0.07(1.7) 

20095-01-18-00 

5.9 (1.6) 

B 

50613.800-50614.155 

0.36 (8.5) 

20094-01 -02- [04,040] c 

21.6 (6.0) 

C 

51603.985-51604.316 

0.33 (8.0) 

40059-01-01-[00,02] c 

13.1 (3.6) 

D 

52615.803-52618.276 

2.47 (59.4) 

70020-01-[02-(01:02),03-01,03-04,04-(00:04)] 

37.9(10.5) 

E 

53013.822-53017.845 

4.02 (96.6) 

70020-03-04-[00:20] 

34.6 (9.6) 

F 

53163.821-53165.483 

1.66 (39.9) 

80027-02-[02-02,02-06,03-(00:02),03-06] 

31.0 (8.6) 

G 

51831.981-51832.459 

0.48(11.5) 

50136-01-04-[00:06] 

18.7 (5.2) 

H 

53168.550-53168.789 

0.24 (5.7) 

80027-02-04-01 

14.1 (3.9) 

I 

52787.744-52787.894 (R) 

0.15 (3.6) 

80114-04-01-[02:04] 

4.4 (1.2) 


52951.992-52952.229 

0.24 (5.7) 

80027-01-01-[02:03] 

9.4 (2.6) 

J 

53003.507-53003.590 

0.08 (2.0) 

70020-03-01-01 

5.6 (1.6) 

K 

53271.381-53271.787 

0.41 (9.8) 

90025-01-02-[02,24,25,27] 

7.9 (2.2) 

L 

55343.361-55350.716 

7.36 (176.5) 

95422-01-[03-(03:04),04-(00:13)] 

49.4(13.7) 


Notes. Subsets colored in red/green in Figure[5]are denoted by (R)/(G) in the Time Interval column. In the ObsIDs column a colon denotes a range. 
a The interval length is shown in units of both days and hours. 
b The total exposure time is shown in units of both ks and hours. 

c The early parts of observations 20094-01-02-040 and 40059-01-01-00 were omitted due to secular motion. 



Soft color Intensity (counts s 1 PCU ’) 


Figure 10. CD/HID showing two instances of Cir X-l undergoing an atoll 
transition from the soft to hard state. Black data points are from a day 
period in 2010 May/June (also shown in panel L in Figure [9} and red data 
points are from a ~7 day period in 2010 August. The data were binned with a 
minimum of 128,000 counts per data point. In some cases entire observations 
did not have that many counts and were combined with other ones close in 
time. 


higher intensities (panels A-C). Going from panel D to A in 
the Cir X-l sequence, the overall color values of tracks de¬ 
crease and the color range spanned drops sharply as the in¬ 
tensity swings become larger, resulting in very stretched and 
flat tracks in the HID. In these “extreme” Cyg-like Z tracks 
in panels A and B Cir X-l exhibits much lower color values 
than ever observed for any of the other three sources analyzed 
in this paper, and it is striking how small the color variations 
along the tracks are despite the large changes in intensity. 


3.4. GX13+1 

GX 13+1 is a bright X-ray binary whose classification as a 
Z or atoll source is ambiguous^ although it has usually been 
labeled an atoll source. [Hasinger & van der Klis ( |1989| ) classi¬ 
fied GX 13+1 as a bright atoll source based on CDs and power 
spectra from EXOSAT observations. In RXTE observations 
from 1996, Homan et al. 1 ( [ 19981 ) discovered a 57-69 Hz QPO, 
which showed similar behavior to the horizontal-branc h os¬ 
cillation (HBO) seen in Z sources. Schnerr et al. p003) per¬ 


formed a combined CD, HID, and power-spectral analysis of 
a large number of RXTE observations made in 1998. They 
found that the source traced out, on a timescale of hours, a 
curved two-branched track in the CD, which showed strong 
secular motion on a timescale of ~1 week. The shape of the 
track was similar to the lower part (IS, LB, UB) of an atoll 
track or the NB/FB part of a Z track; the location of the ver¬ 
tex between the two branches was seen to approximately fol¬ 
low a straight line in the CD. They also found that the source 
showed peculiar CD/HID and rapid-variability behavior com¬ 
pared to most other Z or atoll s ources , but overall they favored 
an atoll classification. Homan et al.| (j2004| analyzed two si¬ 
multaneous A ALL/radio observations of GX 13+1 performed 
in 1999. Based on the results of spectral fits, rapid-variability 
properties, behavior in the radio band, and the scarcity of 
type I X-ray bursts observed from the source since its dis¬ 
covery, they concluded that the properties of GX 13+1 were 
more similar to Z sources than atolls. 

3.4.1. Analysis 

In our analysis we considered all 92 individual ObsIDs for 
pointed observations of GX 13+1. These span a period of 
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Figure 11. CD and HID representing the entire RXTE PC A data set of 
GX 13+1. 


~14 years. In Figure 11 we show a CD/HID based on all 
the GX 13+1 data; strong secular motion is apparent. How¬ 
ever, we note that GX 13+1 overall shows the smallest range 
in secular evolution among the four sources studied here (e.g., 
as quantified by the range in soft or hard color over which the 
lower vertex is observed), and the other three sources show 
secular motion on shorter timescales than does GX 13+1. As 
for Cyg X-2 and Cir X-l we organized the data into subsets; 
this resulted in ~50 such sets, which could in a few cases span 
intervals as long as several days with little or no visible secu¬ 
lar motion. We constructed a sequence of six CD/HID tracks 
(shown in Figure |T2j), which illustrates the secular evolution 
of the source. Similar to our analysis of Cyg X-2 we com¬ 
bined subsets widely separated in time in three of these tracks 
to create more complete tracks than otherwise available and 
thereby give a fuller depiction of the overall secular behavior 
of the source. For these three tracks we also indicate in red 
a representative subset obtained in a short time interval. The 
number of individual subsets used for each track ranges from 
1 to 4; these are listed in Table[5]along with the corresponding 
ObsIDs. Similar to our Cyg X-2 and Cir X-l analysis we or¬ 
der the tracks based on the location of the lower vertex in the 
CD, starting at the lowest soft and hard color. (See discussion 
of this vertex in the following paragraph.) The data used in 
the six tracks together constitute ^42% of the total exposure 
time of the GX 13+1 data set. We note that in a few of the 
subsets used we do see indications of small secular shifts— 
especially in panels B and F in Figure [12] These shifts are 
apparent mostly in the HID (and SID), rather than the CD. 
However, the overall conclusions we draw from these tracks 
and the sequence as a whole are not affected by this. 

As can be seen from Figure [12] the tracks in the CD 
mostly have a two-branched form. The vertex between these 
branchesfollows rather closely a straight line (as observed by 
|Schnerr et al.||2003|), which we show in the figure. We clas¬ 
sify GX 13+1 as a Z source based on its overall secular evo¬ 
lution in the CD/HID, the timescales on which it traces out 
its CD/HID tracks, and its rapid-variability properties. (We 
discuss this further in Section [46] ) We identify the vertex in 
the CD as the lower Z track vertex and the branches to the 
left/right of the vertex line in the CD as the normal/flaring 
branch. In addition, the tracks in panels D and E include 
subsets that we identify as the HB (plus upturn). In the CD 
these segments look very similar to excursions onto the HB 
and HB upturn in the Cyg-like Z tracks of XTE J1701-462, 
Cyg X-2, and Cir X-l. In the HID these segments also stand 
out and are located above and to the left of the rest of the 
track. We note that there is considerable ambiguity regarding 


where within panels C-E to place these presumed HB seg¬ 
ments. What we show in Figure [l2|is therefore our best guess 
but only one of a few possible ways of incorporating them into 
the sequence. What does seem clear from inspection of the en¬ 
tire GX 13+1 data set is that these segments cannot plausibly 
be combined with other track segments—or accommodated in 
our sequence—as NB segments similar to the other ones we 
see. (We discuss these tracks further in Section]?]). 

The tracks in the HID look quite different from those of 
the other three sources. As we discuss below, different track 
branches can be identified in the HIDs; however, the fact 
that these branches are in general very broad and rather ir¬ 
regularly shaped makes them less useful for judging in what 
cases observations can be appropriately combined to form 
CD/HID tracks without significant secular shifts. The SID— 
although also displaying broad and somewhat irregular track 
branches—turned out to be more useful in this respect. To bet¬ 
ter illustrate the behavior of the source in the HID we show in 
Figure |T3] a color-coded version of the track in panel C. The 
locations of the NB/FB vertex in the HID are in general less 
well defined than in the CD and they do not seem to follow a 
straight line. 

3.4.2. Secular Evolution 

GX 13+1 shows behavior that is in some ways quite dif¬ 
ferent from that of the other three sources, especially in the 
HID. A more detailed description of the secular evolution of 
this source is therefore warranted. We note at the outset that, 
as Figure [TT| indicates, the largest portion of the total expo¬ 
sure time for the source was spent on tracks near the middle 
of the sequence in Figure [12] (especially panel D), and it is 
therefore natural that the tracks in these panels would be the 
most complete ones. Conversely, the amount of data available 
at the lowest (panels A and B) and highest (panel F) vertex 
locations (all of which is shown in Figure [12] ) is very small. 
These tracks therefore seem most likely to suffer from incom¬ 
pleteness. 

We first focus on the evolution of the tracks in the CD. As 
the tracks move up the vertex line the FB gradually rotates 
counterclockwise and becomes longer and straighter (more 
Sco-like). Starting at the other end of the sequence, the NB 
gradually becomes longer as the tracks move down the vertex 
line from panel F to D. The NB seems to shorten again going 
from panel C to A, although we suspect that this may be due 
to incompleteness in panels A, B, and possibly C. These three 
tracks may also be missing HB and HB upturn portions. 

The evolution in the HID is less obvious. At the lower ver¬ 
tex locations the NB is strongly tilted to the right and for a 
given hard-color value spans a large intensity range, espe¬ 
cially in its upper part. The NB gradually becomes more 
compact going to higher vertex locations (although it is al¬ 
ways broad compared to the NBs observed in the other three 
sources). The behavior of the FB in the HID is even harder 
to discern than that of the NB. The morphology of the FB in 
panel C can more easily be seen in the color-coded version in 
Figure [13] The upward-pointing part of the FB is essentially 
verticaliffut very short) in panel C, and is at the lowest inten¬ 
sities. Going counterclockwise along this track in the CD, the 
highest intensities are near the top of the NB and then the in¬ 
tensity becomes gradually lower moving along the track until 
the FB reaches its lowest point (in hard color), after which the 
intensity becomes approximately constant. (The same behav¬ 
ior is seen in panel B.) In panels D and E, the flat part of the 
FB in the CD forms a broad patch at the bottom of the track 
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Figure 12. A sequence of CDs/HIDs for GX 13+1 illustrating the secular evolution of the source. Table [5] indicates the data used in each panel. In panels C, 
D, and E data from observations widely separated in time were combined; in each of these three panels a subset of the data—obtained over a period of at most 
~1.8 days—is shown in red. The dashed line in the CDs shows the approximate path followed by the NB/FB vertex. 


Table 5 

Time Intervals and Observations Used to Create GX 13+1 Tracks 




Interval Length 


Exp. Time 

Panel 

Time Interval (MJD) 

days (hr) a 

ObsIDs 

ks (hr) b 

A 

51278.745-51278.971 

0.23 (5.4) 

40023-03-02-03 

10.3 (2.9) 

B 

51276.746-51278.117 

1.37 (32.9) 

40023-03-[01-00,01-02,02-(00:02),02-04] 

26.0 (7.2) 

C 

50990.407-50990.711 

0.30 (7.3) 

30051-01-09-[00:01] 

9.8 (2.7) 


53767.006-53767.601 (R) 

0.60 (14.3) 

91007-08-02- [00,000] 

29.5 (8.2) 

D 

50981.074-50982.844 (R) 

1.77 (42.5) 

30050-01-01-[04:08,050,080] 

54.4(15.1) 


50984.407-50984.576 

0.17 (4.1) 

30050-01-02-03 

8.1 (2.3) 


51390.995-51394.488 

3.49 (83.8) 

40022-01-[01-(00:01),02-00,02-000] 

39.0(10.8) 


51447.322-51447.428 

0.11 (2.5) 

40023-03-04-00 

6.8 (1.9) 

E 

51007.409-51007.560 

0.15(3.6) 

30051-01-12-01 

8.0 (2.2) 


53409.846-53409.949 (R) 

0.10(2.5) 

90173-01-01-00 

6.1 (1.7) 


54740.627-54743.339 

2.71 (65.1) 

93046-08-[02-00,03-(00:01)] 

7.8 (2.2) 


55409.062-55413.805 

4.74(113.8) 

95338-01-[01-(00:07),03-00,03-05] 

24.3 (6.8) 

F 

50997.607-51003.646 

6.04 (144.9) 

30051-01-[ll-(00:03),12-00] 

17.5 (4.8) 


Notes. Subsets colored in red in FigureQ2]are denoted by (R) in the Time Interval column. In the ObsIDs column a colon denotes a range. 
a The interval length is shown in units of both days and hours. 
b The total exposure time is shown in units of both ks and hours. 


in the HID. The rising part of the FB extends from the left 
side of this patch, tilted slightly to the right from vertical, and 
overlaps with the broad NB (see, e.g., the red-colored data in 
panel E). 

3.4.3. Comparison with XTE J1701^162, Cyg X-2, and Cir X-l 

The tracks in the CD of GX 13+1 show many similari¬ 
ties to those of XTE J1701-462, Cyg X-2, and Cir X-l. The 


evolution of the FB in the CD is similar to that seen for 
XTE J1701-462 in selections A-H—with the FB rotating 
counterclockwise and growing longer (and more Sco-like)— 
and likewise similar to that of Cir X-1. There are strong sim¬ 
ilarities between the shape of the FB of GX 13+1 in panels 
D-F and, e.g., the FB of XTE J1701-462 in selections E- 
H, Cyg X-2 in panel J, and Cir X-l in panels H and I. As in 
the other three sources, a lower vertex in the CD of GX 13+1 
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Figure 13. Color-coded version of the track in panel C in Figure [12] illus¬ 
trating the portions of the HID track corresponding to several segments along 
the track in the CD. 


moves along a straight line with positive slope. The NB is also 
similar to the ones of the other three sources (being mostly 
straight and oriented up and to the right), and the possible HB 
and HB upturn segments in GX 13+1 have a similar appear¬ 
ance to those of the other sources as well. 

Although less obvious at first sight, there are also similari¬ 
ties in the HID between GX 13+1 and the other three sources. 
The NB is oriented up and to the right, becoming gradu¬ 
ally shorter and closer to vertical as the tracks move up in 
hard color—i.e., the intensity swings along the NB become 
smaller—and at the same time the maximum intensity ob¬ 
served on the NB gradually decreases (although for this to 
extend to panels A and B we need to assume they have incom¬ 
plete NBs). This is also the case for the other three sources. 
What we interpret as traversals onto the HB and HB upturn 
in GX 13+1 are manifested as movement toward lower inten¬ 
sities and higher hard color; this is in general the case for Z 
sources, including the other three discussed in this paper. In 
the earlier panels in the GX 13+1 sequence the initial part of 
the FB is toward lower intensities (i.e., dipping FB behavior), 
as is seen in the other three sources. The GX 13+1 FB ro¬ 
tates slightly clockwise in the HID (similar to Cyg X-2) as the 
source moves up the vertex line in the CD, evolving into an FB 
that shows an intensity increase rather than decrease. A no¬ 
table difference between GX 13+1 and the other three sources 
is the movement of the lower vertex in the HID. The loca¬ 
tion of the lower vertex in the HID of GX 13+1 is in general 
rather poorly defined, and does not seem to follow a straight 
line. However, it seems clear that the lower vertices in pan¬ 
els A and B in Figure [12] are at a significantly lower intensity 
(^550-600 counts s _ ^PCU _1 ) than the ones in panels C-F 
(^800-900 counts s -1 PCU -1 ), in contrast to the other three 
sources, where in general the intensity at the lower vertex de¬ 
creases as it moves up the vertex line in the CD. 

As is clear from the above GX 13+1 shows similarities to 
both Cyg-like and Sco-like Z sources. Overall, the earlier pan¬ 
els in the track sequence are more Cyg-like and the later ones 
more Sco-like, although the distinction between the two is less 
clear than for the other three sources. 

4. DISCUSSION 
4.1. Secular Evolution 

The main goal of this paper is to study secular evolution 
in the CDs and HIDs of NS-LMXBs, using three sources 
historically known to show substantial changes in the shape 
and location of their CD/HID tracks: Cyg X-2, Cir X-l, and 
GX 13+1. In particular, we want to determine to what ex¬ 


tent the secular evolution of these three sources is similar to 
that seen in XTE J1701-462, the first source found to evolve 
through all NS-LMXB subclasses. 

In Figure [14] we provide an overview of the secular evo¬ 
lution of theCD/HID tracks of XTE J1701-462, Cyg X-2, 
Cir X-l, and GX 13+1. In each of the four panels we show all 
tracks from the sequences in Figures|2][5][9] and[l2]in a single 
CD/HID. (For XTE J1701-462 we show only the red-colored 
subsets from Figured where applicable.) Strong secular evo¬ 
lution is found in all four sources, consistent with reports in 
the literature. While for XTE J1701-462 all data come from 
the dense monitoring of a single 19 month outburst, the data 
for the other three sources were collected sporadically over a 
time span of ~ 14— 15 years. As a result, the secular evolu¬ 
tion in these sources could not be followed in “real time,” as 
was possible for XTE J1701-462, but had to be reconstructed 
from multiple isolated (partial) tracks. To some extent this 
limited our ability to compare the secular evolution between 
sources, in particular for GX 13+1. Despite these issues, it is 
obvious from Figure [14] that the secular evolution of the four 
sources has many common characteristics. 

In the following we summarize our key findings, focusing 
on the similarities between the systems. 

1. As part of their secular evolution, we see clear and con¬ 
tinuous transitions between different NS-LMXB sub¬ 
classes in all four sources. While some of the behav¬ 
ior we report has been described in previous works, 
our work for the first time unambigously links strong 
secular evolution to transitions between various NS- 
LMXB subclasses in Cyg X-2, Cir X-l, and GX 13+1. 
Of the four sources XTE J1701-462 and Cir X-l have 
shown the largest range in behavior; they have dis¬ 
played Cyg-like and Sco-like Z tracks, atoll soft and 
hard states, and have both at some point entered quies¬ 
cence or near-quiescence (^10 35 erg s -1 in the case of 
Cir X-l; |Sell et al 12010) . Cyg X-2 and GX 13+1 have 
only shown Cyg-like and Sco-like Z tracks, with the 
secular changes in GX 13+1 being the most moderate 
of the four sources. 

2. Cyg-like Z source behavior (with large variations in in¬ 
tensity along the NB and HB, and a “dipping” FB that 
is directed toward lower intensities) is observed at the 
highest overall intensities, Sco-like Z source behavior 
(with small intensity variations along the NB and, if 
present, the HB, and large increases in intensity along 
the FB) at lower intensities, and (in XTE J1701-462 
and Cir X-l) atoll behavior at the lowest overall inten¬ 
sities. 

3. Although possible incompleteness in some of the 
CD/HID tracks clouds the picture somewhat, the order 
in which the Z source branches evolve and disappear 
in Cyg X-2 and Cir X-l seems to be consistent with 
what is seen in XTE J1701-462. As the overall inten¬ 
sity decreases, the HB shortens and disappears first, fol¬ 
lowed by the NB, and finally the FB. The situation is 
less clear for GX 13+1; the earlier tracks in Figure [12] 
would need to be missing an HB (and possibly part or 
the NB) for the behavior of GX 13+1 to be consistent 
with XTE J1701-462 as well. 

4. As the shapes of the Z tracks and branches grad¬ 
ually change, the vertices of the tracks shift along 
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Figure 14. CDs/HIDs illustrating the overall secular evolution of XTE J1701-462, Cyg X-2, Cir X-l, and GX 13+1. In the panels we show (from top to bottom) 
all the tracks from Figures [2][5][9] and[l2] respectively. Different colors are used to distinguish between the individual tracks in each panel. 
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well-defined lines in the CD, moving toward higher 
color values as the overall intensity decreases. (We 
note, however, that the overall intensity evolution is 
less clear-cut for GX 13+1, and that CirX-1 shows 
some deviations from the general intensity trend.) For 
XTE J1701-462 and Cyg X-2 this also applies to the 
HID, and to a large extent for Cir X-l as well. 

5. As the overall intensity of the Z tracks increases and 
they become spectrally softer, their dynamic range (i.e., 
ratio of maximum to minimum intensity) increases as 
well. Whether this trend also extends to the tracks at 
the lowest color values in GX 13+1 is unclear due to 
possible incompleteness issues. 

Despite the fact that differences exist in their detailed be¬ 
havior, based on the above findings we conclude that the sec¬ 
ular evolution in the four NS-LMXBs that we studied largely 
follows similar patterns. This strongly suggests that the se¬ 
quences of NS-LMXB subclasses reported in this work are 
representative of the class of NS-LMXB as a whole. 


4.2. Evolution of the Flaring Branch 

Of the three main Z source branches, the FB undergoes the 
most dramatic changes in its shape and orientation. In all four 
sources it rotates in both the CD and HID and evolves from a 
branch that shows a strong decrease in intensity with respect 
to the NB/FB vertex (in the Cyg-like Z tracks), to one that 
shows a strong increase in intensity with respect to that vertex 
(in the Sco-like Z tracks), although the amplitude of these in¬ 
tensity changes is significantly smaller in GX 13+1 than in the 
other three sources. From the dense “coverage” of the secular 
evolution of XTE J1701-462 and Cyg X-2 it is clear that this 
FB evolution from dipping to flaring is a very gradual process, 
with the morphology changing smoothly throughout the track 
sequences. This suggests that the dipping Cyg-like FB and the 
Sco-like FB are related phenomena, and it seems unlikely that 
the overall observed FB behavior can be explained by two pre¬ 
sumably unconnected mechanisms: absorption by the outer 
disk, which was proposed as an explanationfor the Cyg X-2 
dipping behavior by |Ba!ucinska-Church et al. ( 2011| ), and nu¬ 
clear burning on the neutron star (combined with increases in 
the mass accretion rate in the Sco-like Z s ources), pr oposed 
as an explanation_for the flaring FB by [ Church et al. ( |2012| . 
(See also |BaIucinska-Church et al.|2012[ ) 


4.3. Other Z Source Transients 

We note that two additional Z source transients discovered 
after we began work on this paper, IGR J17480-2446 (also 
known as Terzan 5 X-2) and MAXI J0556-332, have shown 
similar behavior to the sources studied here. RXTE observa¬ 
tions of IGR J17480-2446 revealed clear Cyg-like Z source 
behavior at the peak of its ^2.5 month outburst and atoll be¬ 
havior at lower intensities; due to the rapid source evolution 
and sparsity of the data Sco-like Z behavior could not be iden¬ 
tified unambiguously {A ltamira no et al.||2010| [Chakraborty 
|et al.|20lf| |Altamirano et al.|2012| D. Altamirano et al. 2015, 
in preparation). 

MAXI J0556-332 showed both Cyg-like and Sco-like Z 
behavior during its 16 month outbur st ( Homan et al.||20ir 


Sugi zaki et al. 2013| |Homan et al.||2014| J. Homan et al. 
2015, in preparation). Figure 2 in Homan et al. | ( |20 1 4] ) shows 


a combined HID using all the RXTE data for the source 


(panel a) as well as four time-selected HID tracks illustrat¬ 
ing the secular evolution of the source (panels b-e). The 
total HID looks similar to that of Cyg X-2 (Figure [ 3 ]). The 
higher-intensity MAXI J0556-332 tracks (panels b and c) are 
Cyg-like Z tracks; the track in panel b shows an extended 
dipping FB and resembles the highest-intensity tracks of 
XTE J1701-462, Cyg X-2, and Cir X-l. The lower-intensity 
tracks of MAXI J0556-332 (panels d and e) are Sco-like Z 
tracks, resembling the tracks in panels F and G-H, respec¬ 
tively, in the XTE J1701-462 sequence, and K+L and M in the 
Cyg X-2 sequence. As the overall intensity changes and the 
MAXI J0556-332 tracks gradually evolve in shape, both the 
upper and lower vertices follow straight lines in the HID, as 
observed for XTE J1701-462 and Cyg X-2. The RXTE cov¬ 
erage of MAXI J0556-332 ended before the outburst came to 
an end, while the source was still showing Sco-like Z behav¬ 
ior. It may subsequently have passed through the atoll soft 
and hard states on its way back to quiescence. 

Finally, we note that the atoll transient XTE J1806-246 
likely entered the Sco-like Z source regime briefly at the peak 
of its outburst in 1998 ( [Wijnands & van der Klis|1999| ). Ob¬ 
servations during the outburst peak show the source tracing 
out a curved (partial) track resembling the lower NB+FB part 
of a Sco-like Z track. As one would expect for a transition 
from the atoll to Z regime, this track is at lower color values 
(and higher intensities) with respect to the observations in the 
atoll soft state taking place before and after the outburst peak. 
While the source traced out this track, a 7-14 Hz QPO was 
detected, whose frequency increases as a function of position 
along the track (going from the NB to FB), consistent with the 
behavior of the normal/flaring-branch oscillation in Z sources. 

The three sources discussed in this section, along with the 
four sources analyzed in this paper, constitute the entire sam¬ 
ple of NS-LMXBs known or believed to have shown transi¬ 
tions between at least two of the three main NS-LMXB sub¬ 
classes (Cyg-like Z, Sco-like Z, atoll). The fact that these 
three additional transients show behavior seemingly largely 
consistent with the other four sources (e.g., the relative lumi¬ 
nosity dependence of the subc lasses) strengthens our conclu¬ 
sion at the end of Section ItTI 


4.4. Role of the Mass Accretion Rate 

The spectral fitting results of |Lin et al.| ( 20091)] ) indicate 
that the secular evolution of XTE J1701-462 was driven by 
changes in the mass accretion rate—they find that the M in¬ 
ferred from the accretion disk component in the spectra at 
the upper and lower vertices decreases monotonically when 
moving left along the two vertex lines in the HID. (See also 
discussion in |H10[ ) Based on the similarities in secular be¬ 
havior between XTE J1701-462 and the three sources stud¬ 
ied in this paper we conjecture that the secular evolution of 
Cyg X-2, CirX-1, and GX 13+1 is also driven by changes 
in the mass accretion rate. The case is clearest for Cyg X-2, 
whose behavior is most regular of the three sources and clos¬ 
est to that of XTE J1701-462. The fact that the overall in¬ 
tensity of the Cyg X-2 tracks decreases monotonically as they 
gradually evolve in shape and move along the vertex lines to 
higher color values makes it plausible that the mass accretion 
rate in Cyg X-2 decreases monotonically along the sequence, 
as seems to be the case for XTE J1701-462. 

Although there are some important differences between the 
evolution of Cir X-l and XTE J1701-462, the global proper¬ 
ties are still very similar. The overall large intensity decrease 
from the Cyg-like tracks in the earlier panels to the Sco-like 
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(and finally atoll) tracks in the later panels suggests a decreas¬ 
ing mass accretion rate going from panel A to L. However, the 
weaker correlation between track shape and location in the 
HID compared to the other sources—in particular manifested 
in the jumping back and forth in the HID seen during part of 
the sequence (panels F-I)—suggests that the mass accretion 
rate progression m ay n ot be strictly monotonic along the se¬ 
quence. In Section H3] we suggest that these “discrepancies” 
between the track shape and location in the HID may be the 
result of rapid changes in the mass accretion rate. 

Despite the differences in the behavior of GX 13+1 com¬ 
pared to the other three sources we find it most likely that 
changes in the mass accretion rate are also responsible for 
the secular evolution of GX 13+1 and that, as in the other 
sources, M decreases as the source moves up the vertex line 
in the CD—i.e., going from panel A to F in Figure fl2| As 
the tracks become overall spectrally harder the maximum in¬ 
tensity of a given track decreases, and the intensity swings 
along the track become smaller, as observed for the other three 
sources, although incompleteness in the tracks in panels A and 
B needs to be invoked to explain their deviations from this 
trend. Such incompleteness can, however, not explain the fact 
that the lower vertex in the HID is at lower intensities in those 
two tracks than in tracks that we suggest exhibit lower mass 
accretion rates, opposite to the behavior of XTE J1701-462. 

Our results support the idea that the mass accretion rate 
is responsible for changes between subclasses in individual 
sources, with Sco-like Z behavior appearing at lower accre¬ 
tion rates and overall luminosities than Cyg-like Z behavior 
(and atoll behavior appearing at accretion rates that are lower 
still). However, in this picture other parameters, such as the 
neutron star spin, mass, and magnetic field, may still affect 
the luminosity scale of the subclass sequence—i.e., the range 
in luminosity (and mass accretion rate) over which a given 
s ubclass appears may differ between sources. 

|Lin et al.j ( [2009a|b| ) showed that for its best-estimate dis¬ 
tance of 8.8 kpc, XTEJ1701-462 reached super-Eddington 
luminosities during its Cyg-like stage. Even for a distance 
that is 30% smaller this would have been true. The super- 
Eddington nature of the Cyg-like tracks is further supported 
by observations of two Sco-like Z sources with accurately 
know n distances, ScoX-1 (2.8± 0.3 kpc ; Bradshawet^a 


1999 ) and LMCX-2 (50±2 kpc; |Pietrzyns ki et al. 2013 
Spectral fits indicate that those sourc es are either around Ec 


dington or mildly supe r-Eddington (Bradshaw et al.| 2003 


Agrawal & Misra 2009). Given that Cyg-like behavior is ob 
served at significantly higher overall intensities than Sco-like 
behavior in XTE J1701-462, Cyg X-2, and Cir X-l, this sug¬ 
gests that Cyg-like Z tracks may in general exhibit luminosi¬ 
ties well above the Eddington limit. 

As the sources analyzed in this paper evolve from Sco-like 
to Cyg-like Z behavior, it is interesting to see a large increase 
in the dynamic range (i.e., ratio of maximum to minimum 
intensity) of the Z tracks. This is most clearly seen in the 
tracks of XTE J1701-462, Cyg X-2, and Cir X-l, although 
this also holds for GX 13+1, apart from the possibly very in¬ 
complete track in panel A (Figure 12). The most extreme in¬ 
tensity swings are seen in Cir X-l, where we observe changes 
in the intensity (going between the tip of the dipping FB and 
the upper vertex) by factors up to 8 in periods as short as 
20 minutes in the highest-intensity tracks. In these extreme 
Cyg-like Z tracks, Cir X-l seems likely to have reached sig¬ 
nificantly higher luminosities (perhaps by a factor of ^3) and 
exhibited significantly higher mass accretion rates than ob¬ 


served in XTE J1701-462 at the peak of its outburst. Given 
that the large Cyg-like intensity swings observed in Cir X-l 
(as well as XTE J1701-462 and Cyg X-2) occur in inferred 
luminosity ranges that are near- or super-Eddington, it seems 
likely that radiation pressure effects play an important role. 
While it is unlikely that the mass accretion rate toward the in¬ 
ner parts of the accretion flow varies this strongly on such a 
short timescale, strong radiation pressure may result in sub¬ 
stantial changes in the geometry of the inner accretion flow, 
with a (large) fraction of the mass inflow possibly being redi¬ 
rected into wind or jet outflows, resulting in large intensity 
swings. 

Given the near-/super-Eddington nature of the Z sources, 
it is also interesting to compare them with the transient ul- 
traluminous X-ray source M82 X-2. This source was re¬ 


cently identified as an accreting neutron star (Bachetti et al. 
2014) and has been s een to reach a ma ximum luminosity 
of ~1.8xl0 40 erg s -1 (Feng et al. 2010), which is a factor 
of ^30 higher than the peak luminosity of XTE J1701-462 
(Lin et aL| |2009b) . However, the short-term light curves of 


M82 X-2 presented in |Kong et al.| ( |2007'] ) do not show any of 
the strong luminosity swings that characterize the Cyg-like 
Z sources. A possible explanation for this is that due to the 
much higher neutron star magnetic field in M82 X-2 (as in¬ 
ferred from the strong pulsations in the light curve; Bachetti 
et al. 2014| ) the accretion flow geometry in that source isvery 
different from that in the Cyg-like Z sources, perhaps prevent¬ 
ing the strong feedback mechanisms that we think may cause 
the large luminosity swings. 

Three of the sources studied in this paper (XTE J1701-462, 
Cyg X-2, and GX 13+1) were also part of a large-sampl e 
study of NS-LMXB behavior by Munoz-Darias et al.|([2014|. 
These authors find that the Cyg-like and Sco-like Z sources 
occupy a region of the rms-luminosity diagram that partly 
overlaps with that of the luminous black hole X-ray binary 
GRS 1915+105. We note, however, that overall the Z sources 
and GRS 1915+105 behave very differently (e.g., Belloni| 
et al.||2000| ), suggesting that at near/super-Eddington lumi- 
nosities the observed properties of NS-LMXBs and black hole 
X-ray binaries are quite different, perhaps (in part) due to the 
increased importance of radiation feedback from the neutron 
star surface. 

4.5. Nonstationary Accretion in Cir X-l 


As discussed in Sections |3.3| and |4.4| the secular behavior 
of Cir X-l is not as regular as that of XTE J1701-462 and 
Cyg X-2, as for example manifested in deviations of the track 
vertices from a straight line path (especially in the HID). We 
suspect that this may be the result of the highly nonstationary 
accretion in Cir X-l, thought to arise from the binary’s ec¬ 
centric orbit. From transient LMXBs it is known that highly 
nonstationary accretion can lead to substantial hyste resis (e.g., 
Maccarone & Coppi|2003[ [Homan & Belloni|2005] ), with cer- 
tain source states being observed over a much larger luminos¬ 
ity range than would be the case for slowly changing accre¬ 
tion rates. [Yu & Yan| ( |2009| showed that these effects become 
stronger as the changes in mass accretion rate become faster. 
We speculate that, analogously, highly nonstationary accre¬ 
tion could be a possible explanation for certain Z track shapes 
being observed over a range of luminosities, rather than at a 
single location along the vertex lines. While all four sources 
in our sample undergo nonstationary accretion, which gives 
rise to their secular evolution, the rate at which the mass ac¬ 
cretion rate varies with time is likely considerably higher in 
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CirX-1 than in the other three sources. This is suggested 
by the fastest timescales on which we observe secular evo¬ 
lution in Cir X-l, which are shorter than in the other sources. 
The fastest secular evolution seen for XTE J1701-462 took 
place during the earliest stages of the outburst (panels A-C). 
Inspection of the SID indicates that the lower vertex of the 
HID track in panel A is at ^1900 counts s -1 PCU -1 rather 
than the ^1600 counts s -1 PCU -1 implied by the vertex line 
shown; we see indications for a similar, but smaller, shift in 
panel B. We also note that MAXI J0556-332 showed strong 
deviations from a straight vertex line (which it otherwise fol¬ 
lowed) during the rapid rise of its outburst (J. Homan et al. 
2015, in preparation). This lends support to our hypothesis 
that rapid changes in the mass accretion rate underlie the ob¬ 
served deviations in Cir X-l. 

An additional factor that may play a role in the somewhat 
irregular behavior of Cir X-l is the fact that the neutron star 
spin, the companion spin, and the orbital axis are likely mutu- 
ally misaligned in this system (Brandt & Podsiadlowski 1995 


|Heinz et al.|2013| ). This can lead to precession of the neutron 
star spin axis, which may in turn affect the rotation of the in¬ 
ner accretion flow. Changes in the projected area of the inner 
accretion flow might then result in similar track shapes being 
observed at different X-ray fluxes. 


4.6. The Z Source Classification of GX13+1 

Overall, we find that the shape and evolution of the CD/HID 
tracks of GX 13+1 show much stronger similarities to the 
other three sources analyzed in this paper than to atoll sources. 
Most generally, just the fact that we see gradual evolution in 
the shape of large-scale tracks traced out by the source as the 
locations of the tracks shift over a considerable range in the 
CD and HID points to a Z source, since such secular evolution 
in the CD/HID has never been observed in atoll sources—only 
in Z sources. Moreover, the specifics of this evolution show 
many similarities to those of X TE J 1701-462, Cyg X-2, and 
Cir X-l, as discussed in Section [34] Furthermore, our identi¬ 
fication of track segments in the GX 13+1 sequence (panels D 
and E in Figure 


in Section 3.4.1 


T2J ) as the HB and HB upturn (see discussion 
_ is supported by the detection of a highly sig¬ 
nificant ~ 22-2$ Hz QPO in the corresponding timing data of 
these segments (J. K. Fridriksson et al. 2015, in preparation). 
These frequencies are considerably lower than the 57-69 Hz 
QPO found by |Homan et al.| ( p~998| ) in parts of the tr acks that 
we now cla ssify as the upper part of the NB (see also Schnerr 
et al.|2003|). If we interpret the QPO as an HBO, the decrease 
in its frequency as the source moves onto the HB from the 
upper NB is consi stent with what h as been seen in other Z 
s ources (see, e.g., Jonker et al.|2000 and Homan et al.|2 002). 

|Munoz-Darias et al.| ( |2014| ) find that in their source sample 
there are systematic differences between the ranges in fast- 
variability level (broadband fractional rms) shown by the dif¬ 
ferent NS-LMXB subclasses. The Cyg-like Z sources reach 
the highest fractional rms values (up to ~10%—15%), fol¬ 
lowed by the Sco-like Z sources (with peak values in the 
^5%-10% range), while bright persistent atoll sources and 
transient atolls in a_ bright soft state are typically observed 


at <5%. Munoz-Darias et al. (2014) find that the rms val 


ues oberved for XTE J1701 7462 in the different subclasses 
are consistent with this. The highest fractional rms in the Z 
sources is observed on the HB, and we find that on the pre¬ 
sumed HB segments GX 13+1 reaches fractional rms values 
as high as ^9%-10% (J. K. Fridriksson et al. 2015, in prepa¬ 


ration). This seems more consistent with Z than bright atoll 
behavior. We also note that the hard color range over which 
these GX 13+1 HB segments are observed (~0.6-0.8) is con¬ 
sistent with HBs observed in both persistent and transient Z 
sources, but clearly too low to be consistent with the atoll hard 
state ( |Fridriksson|201 1| ). 

Considering, in addition to the above, the timescale on 
which the GX 13+1 tracks are traced out (hours to ^1 day 
instead of the days/weeks timescales usually observed in 
atoll sources), the extreme scarcity of observed type J_ X- 
ray bursts from the source ([ Fleischman 1985; Matsuba et al. 


19951 |Galloway et al.||2008|), and the reported radio b right 

- - 


ness/behavior ( Fender & Hendry 2000; Hom an et al.| [2004) 
we confidently classify GX 13+1 as a Z source, one that 
switches between Sco- and Cyg-like Z behavior. However, 
the behavior of the source in the HID remains somewhat puz¬ 
zling when compared with XTE J1701-462, Cyg X-2, and 
CirX-1—especially the broad appearance of the tracks and 
the fact that in the spectrally harder tracks the NB/FB vertex 
is at higher, not lower, intensities than in the softer tracks. 

This behavior may be connected to the fact that GX 13+1 
is probably viewed at a high inclination angle (perhaps 
60°-80°) as indicated by observations of absorption dips in 
this system on several occas ions ( [Diaz Trigo et al. [20 12 |[D t A 7| 
et a l. | 2014l|Iaria et al.|2014| ). It was recently shown for black 


hole X-ray transients that the s hape of their HID tracks de - 
pends on the inclination angle (Munoz-Dari as et aLl|2013| ). 
While it is not fully understood how inclination affects the 
shape of the HID, it is likely that an inclination dependence 
of the HID is present in NS-LMXBs as well, with different 
physical components, such as the accretion disk and bound¬ 
ary layer, being affected differently by changes in the viewing 
angle. 

5. SUMMARY AND CONCLUSIONS 

It was recently shown that during its 2006-2007 outburst, 
the super-Eddington neutron star transient XTE J1701—462 
evolved through all subclasses of low-magnetic-field NS- 
LMXB s as a result of changes in a si ngle physical param¬ 
eter: the mass accretion rate ( |Lin et al.|[2009bl [hTO] ). The 
main goal of this paper is to study secular evolution in the 
CDs and HIDs of NS-LMXBs and, more specifically, to in¬ 
vestigate whether evolution similar to that of XTE J1701-462 
occurs in other sources. To this end we studied the evolu¬ 
tion of the CD/HID tracks of three sources that show compli¬ 
cated secular behavior—Cyg X-2, CirX-1, and GX 13+1— 
using all RXTE PCA data obtained during the lifetime of the 
mission. No comprehensive study of the CD/HID behav¬ 
ior of these sources—using most or all of the RXTE data 
now available—has previously been performed. We cre¬ 
ated sequences of CD/HID tracks—in many cases carefully 
piecing together data obtained throughout the RXTE mission 
to construct more complete individual tracks than otherwise 
available—that demonstrate the secular progression in the 
track morphologies and locations. In the case of Cir X-l, this 
analysis had to be preceded by the filtering out of data affected 
by intrinsic absorption, which otherwise strongly influences 
the appearance of the CD/HID. 

We find that Cyg X-2 and CirX-1 show strong similari¬ 
ties in their secular behavior to XTE J1701-462; in particu¬ 
lar, both sources display gradual evolution between Cyg- and 
Sco-like Z tracks, with the latter occuring at lower intensities. 
Cir X-l shows especially extreme versions of Cyg-like tracks 
at the highest observed intensities, with very large and rapid 
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intensity swings but only small changes in spectral color. At 
the lowest observed intensities we also see clear transitions in 
the CD/HID from the atoll soft to hard state for Cir X-l. The 
very gradual evolution of the FB morphology that is observed 
(especially in XTE J1701-462 and Cyg X-2), from a dipping 
(Cyg-like) FB to a flaring Sco-like one, suggests that these 
two different forms of the FB are related phenomena. 

Although GX 13+1 shows behavior that is in some ways 
peculiar—especially in the HID—we find that the source 
also displays similarities to XTE J1701-462, Cyg X-2, and 
Cir X-l, and overall we conclude that the properties of 
GX 13+1 are strongly indicative of a Z source. 

The fact that Cyg X-2 and Cir X-l show evolution between 
different NS-LMXB subclasses that is similar to the evolu¬ 
tion seen in XTE J1701-462 lends support to the suggestion 
of|H10|that the behavior of XTE J1701-462 is representative 
of the entire class of NS-LMXBs. It also strengthens their 
conclusion that (at least within individual sources) Cyg-like 
Z behavior takes place at higher luminosities and mass accre¬ 
tion rates than Sco-like Z behavior. Based on the similarities 
to XTE J1701-462 we conjecture that the secular evolution 
of Cyg X-2, Cir X-l, and GX 13+1 is due to changes in the 
mass accretion rate in the systems, and overall our results sup¬ 
port the notion that differences in mass accretion rate are the 
primary factor underlying the existence of the various NS- 
LMXB subclasses. We conclude that Cyg X-2 and Cir X-l— 
like XTE J1701-462—probably reach super-Eddington lumi¬ 
nosities, with Cir X-l likely reaching significantly higher lu¬ 
minosities and mass accretion rates than the other sources. 


We thank the referee for constructive comments. This re¬ 
search has made use of data obtained from the High Energy 
Astrophysics Science Archive Research Center (HEASARC), 
provided by NASA’s Goddard Space Flight Center. 
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